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Cosmic string Y-junctions: a comparison between field theoretic and Nambu-Goto 

dynamics 
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We explore the formation of cosmic string Y-junctions when strings of two different types collide, 
which has recently become important since string theory can yield cosmic strings of distinct types. 
Using a model containing two types of local U(l) string and stable composites, we simulate the 
collision of two straight strings and investigate whether the dynamics matches that previously ob- 
tained using the Nambu-Goto action, which is not strictly valid close to the junction. We find that 
the Nambu-Goto action performs only moderately well at predicting when the collision results in 
the formation of a pair of Y-junctions (with a composite string connecting them). However, we find 
that when they do form, the late-time dynamics matches those of the Nambu-Goto approximation 
very closely. We also see little radiative emission from the Y-junction system, which suggests that 
radiative decay due to bridge formation does not appear to be a means via which a cosmological 
network of such string would rapidly lose energy. 
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I. INTRODUCTION 

Cosmic strings [l|, |2[ may not have played the primary 
role in the seeding of cosmic structure, with inflation ap- 
pearing to have had that function, but they may still 
be important cosmological entities. Observations of, for 
example the cosmic microwave background (CMB) ra- 
diation, merely limit the allowed string tension. And 
current datasets do not do so particularly stringently: it 
is required that the string tension is less than about one 
third of that which would see them dominate the temper- 
ature anisotropies in the CMB [3|, 0, 0, @| ■ Indeed they 
may still prove to make the primary contribution to the 
as-of-yet undetected CMB polarization B-mode @, [a Q • 
Future CMB data, galaxy redshift surveys, gravitational 
wave experiments and gravitational lens surveys promise 
to either greatly tighten the existing constraints, or to 
plausibly detect cosmic strings. 

This is particularly important for (super) string/M- 
theory, since it has recently been realized that fundamen- 
tal superstrings need not be limited to microscopic scales. 
Now these fundamental F-strings, along with other string 
theory entities called D-strings, appear able to play the 
role of cosmic strings [10|, LU| • And these cosmic super- 
strings would have particular properties, for example, be- 
cause of the extra dimensions required by string theory, 
or because (p,q) bound states of p F-strings and q D- 
strings can form, with Y-shaped junctions where they 
unzip into two more basic constituents. As a result, the 
detection of cosmic strings would provide an exciting ob- 
servational window upon string theory. 

There has been a great deal of recent work on the topic 
of cosmic superstrings: from studies predicting their for- 
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mation in brane inflation models [13 . 1 1 31 ] right through to 
their dynamics at late times. However, the latter case is 
not completely understood even for the traditional situa- 
tion of gauged U(l) strings, with there being some ques- 
tion marks over the rate at which strings self-intersect 
and chop off small loops 0, M, [JJ, [13 111 El M, H • 
The more complex case of cosmic superstrings is there- 
fore particularly challenging. Various authors have used 
numerical simulations of field theories to represent cosmic 
superstrings on horizon-size scales, including the use of 
linear sigma models [22| and global SU(2)/Z3 strings [23j |. 
as well as more realistic models involving local strings 
[24 [25j . These have largely addressed the question: do 
cosmic superstring networks evolve in the same manner 
as traditional gauge strings, in that their mean energy 
density scales with the total density of the universe?. The 
concern is that the bound states and Y-junctions would 
slow the sub-horizon decay of the strings, resulting in 
them dominating the universe at late times. That would, 
of course, be in clear contradiction with observation but 
fortunatelythe above simulations, as well as analytical 
modeling 26, 27], suggest that superstring networks may 
exhibit scaling. However, this work is very challenging 
and the issue is not completely resolved. 

Given the complex array of string seen in such horizon- 
scale simulations, it is difficult to understand the micro- 
physics involved in the problem and this is essential for a 
reliable understan ding of the results. However Copeland, 
Kibble and Steer [28|, [29|, hereafter CKS, have recently 
used the Nambu-Goto approximation to study the colli- 
sion between two straight strings, and have shed a great 
deal of light upon Y-junction formation. Unfortunately, 
the Nambu-Goto action assumes that the string separa- 
tion and curvature scale are far greater than the string 
width and it is therefore not strictly valid at the site of 
the Y-junction itself. It cannot include, for example, the 
interaction between the strings, which is of course asso- 
ciated with the formation of the bound state. Indeed, 
the attraction between strings in the vicinity of the Y- 



junction may naively be expected to increase the bridge 
growth rate and to allow bridge formation when it is ruled 
out under the Nambu-Goto approximation. Moreover, it 
has been suggested [26| that the energy liberated by the 
formation of stable composites could be released as radi- 
ation and therefore help prevent the network from dom- 
inating the universe, but the Nambu-Goto action does 
not allow for such radiation. 

Hence in this article we investigate the formation and 
dynamics of Y-junctions using 3D field theory simula- 
tions of a model involving two coupled Abelian Higgs 
models, as introduced by one of us in Ref. [30]. Specif- 
ically we study the formation of bound states when two 
straight strings collide and then we compare our results 
with the analytical predictions from CKS. Reference [3l| 
has previously studied such collisions for a single Abelian 
Higgs model in the Type I regime, when bound states and 
Y-junctions can form, and this has been re-visited in or- 
der to test the CKS results by Ref. [32j. Type I strings 
offer a different type of coalesce in that there is a single 
type of flux present, and have been studied recently for 
theories containing supersymmetric flat directions [331 ] - 
Here our model contains two independent U(l) symme- 
tries, modeling the separate F and D string charges, and 
hence our results are of a different nature to those of 
Ref. [32]. Furthermore, our detailed measurements from 
the collision aftermath are the first to quantitatively in- 
vestigate, not merely whether a composite region forms, 
but also its growth rate and precise dynamics. These are 
also important with regard to the understanding of the 
Y-junctions, as well as to the future application of the 
CKS approach. 

In the next section we discuss the dynamics of Y- 
junctions under the Nambu-Goto approximation, before 
discussing the field theoretic model employed here in Sec. 
IIHI We then discuss our simulation method in Sec. IIVI 
and our qualitative results in Sec. fVj Only with those 
results in hand can we discuss our the methods behind 
our detailed measurements from the simulations or the 
results from them, which form Sees. IVII and FVIII respec- 
tively. Finally we then present our interpretation of these 
results and our conclusions. 
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FIG. 1: The intersection of two strings to form a pair of 
Y-junctions, as seen in the Nambu-Goto picture. The initial 
state of the two strings (top) is with them lying parallel to the 
yz-plane and travelling with velocities ±v in the ^-direction, 
while there is a choice of final state (middle and bottom) 
depending upon how the two strings connect to each other. 
In either case, a bridge string links the two Y-junctions and 
the initial strings have kinks travelling along them. We denote 
the middle case as a cosine linkage, and the bottom case as a 
sine linkage. 



II. NAMBU-GOTO DYNAMICS 

The Minkowski space-time dynamics of three strings 
mee ting at a Y-junction was solved analytically by CKS 
[28l . [291 ] under the Nambu-Goto action. Since no at- 
traction between strings is included in this action, their 
approach was to add a Lagrange multiplier in order to 
constrain the three strings to coincide at the junction. 
Initial conditions may then be chosen such that there 
are two strings lying in the y-z plane as shown in Fig. 
QJtop), each making an angle a to the z-axis and travel- 
ling with velocities v and — v in the x-direction. However, 
the Nambu-Goto action cannot illuminate the creation of 
a composite bridge string upon their intersection and it 



must be inserted by hand how the strings connect to each 
other in the final state, with two such possibilities shown 
in Fig. [T] We refer to these here as cosine and sine con- 
nectivities, since the essential difference between them is 
just the choice of placement for the angles a and 9, and 
correspondingly whether cosine or sine terms appear in 
the primary equation below. 

In general there is a third possibility for the connectiv- 
ity of the strings, in which there is no partner exchange 
and a bridge simply grows between the two initial strings 
[281 ] . However, this is forbidden in the field theoretic 
model that we consider here and we will not discuss it 
further. 

Having chosen the connectivity, the next step in the 
CKS method is to solve for the bridge velocity and ori- 



entation. If the two incident strings have energies per 
unit (invariant) length, jj,\ and /U2, that are equal, then 
the symmetry present in the problem greatly simplifies 
matters. Firstly the two Y-junctions must remain in the 
yz-plane at all times. Secondly, the bridge must lie ei- 
ther along the z-axis (cosine link) or along the y-axis (sine 
link), and hence both the bridge speed and orientation 
are trivial 28]. 

Even in an unsymmetric case (/ii 7^/12), the bridge 
must lie still parallel to the yz-plane due to the symmetry 
between the two Y-junctions and its velocity must be 
parallel to the x-axes [29j. For cosine connectivity, as 
shown in the middle pane of Fig [TJ CKS then derive the 
following equation for the bridge velocity u [291 ]: 



= u 4 S 2 sin 2 a 

,2 



(1) 



u- [R 2 (l-v 2 ) + S 2 {v 2 cos 2 
S 2 v 2 cos 2 a, 



sin a) 



where we introduce the notation R = /^/(/ii + ^2) and 
also 5* = (/ii — /i2)/(/ii + ^2)- If Mi — M2, then S is zero 
and the solution is just u — 0, as noted above. For S 7^ 0, 
this equation always yields one positive root for u 2 , while 
the sign of u matches that of S if string 1 initially had 
positive a;- velocity. 

With u having been determined, the angle 9 between 
the bridge and the z-axis may then be found using [29| |: 
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FIG. 2: The Nambu-Goto predictions for when cosine bridges 
can form in cases A, B and C that we will explore using field 
theory simulations, with sine bridges also included for case C. 
Case A has S = and R = 0.84 while case B has 5" = 0.30 
and R = 0.84. Case C has S = 0.25 but in the field theory 
case the value of R depends on the connectivity due to flux 
cancellation for a cosine bridge (Cc) with R — 0.37 but not 
for a sine bridge (Cs), which then has R = 0.90. Case Cc 
therefore has \S\ > R 2 . Note that Y-junctions under both Cc 
and Cs are possible outcomes for a small region of parameter 
space while for cases A and B there is no overlap between sine 
bridges (not shown) and the cosine bridges. 



Of course, if S = u = 0, then 9 is simply zero and, as 
noted above the bridge lies simply on the z-axis for cosine 
connectivity. 

From these values of u and 9, the rate at which the 
length of the bridge grows may be solved for [29]. For 
S t^ it is convenient to use the invariant half-bridge 
length S3 = h/^fu; where the physical bridge length is 2/3 
and 7,j = l/yl — v? is the usual Lorentz factor. This 
yields the neat causality constraint: S3 < 1; since the Y- 
junction, which moves in the yz-plane at speed dl^/dt and 
at speed u in the cc-direction, cannot traverse at super- 
luminal speeds. A simple expression for S3 was noted in 
Ref. 34j and in our notation this is: 



«3 



7 U cos a — Rj cos 9 
7 cos 9 — R~f u cos a 



(3) 



Although there is no explicit dependence upon S here, 
this is realized through 7„ and 9, and can have a large 
impact on the result. However, if S = then 7^ and cos 9 
are both equal to unity and may be therefore omitted, 
leaving a particularly compact expression [281 ]. 

If sine connectivity is chosen instead, then it is simply 
the case that the replacement a — ► (tt/2 — a) and 9 — > 
(tt/2 — 9) should be made, or equivalently for the above 
equation, the cosines are merely changed to sines. 

Of course, the solution only makes physical sense if 
S3 > and Fig. [Jlshows the region of the a-v plane where 



CKS solutions are possible for the three cases that we 
will explore later using field theory simulations. Case A 
has S = and R — 0.84, and Eqn. ([3]) then implies 
that a cosine bridge can only form if a < arccos(7-R). 
That is, for a given v there is a certain a above which 
a CKS solution is not possible and hence it is implied 
that the strings must simply pass through each other. 
Swapping arccos(7-R) for arcsin(7_R) gives the constraint 
for a sine bridge, which for case A is not shown in the 
figure since it is just a repeat of the cosine case but with 
a — ► (n/2 — a). However, if R was less than cos 45° then 
there would be a range of a for which growing solutions 
exist for both connectivities and the Nambu-Goto action 
would not then reveal which solution would be followed. 
This is something we will explore later, but only for a 
case with S 7^ 0. 

Case B has (almost) the same value of R but also has 
S = 0.3. This yields, in fact, a largely similar situation 
with respect to the figure, although with bridge formation 
possible over a slightly larger area. For l^l > R 2 the sit- 
uation becomes quite different and this is possible under 
case C. Considering firstly only cosine connectivity, the 
plotted case Cc highlights that CKS solutions are permit- 
ted for at high v for all a. However, a subtle point is that, 
in our field theory, flux cancellation (see Sec. [V} occurs 
for the cosine connectivity to yield a low R of 0.37 but 
this does not happen for sine connectivity and hence case 



Cs has a larger R of 0.9. Since the condition I.!?! > R 2 
is not met for Cs, then this last case appears simply like 
cases A and B but with a — > (ir/2 — a), although we 
plot it now since it is not simply a repeat of Cc. Impor- 
tantly, there is a region in which both Cc and Cs solutions 
are possible and hence for case C we will simulate initial 
conditions for which the Nambu-Goto dynamics make no 
prediction as to which solution will actually be followed. 
When bridge growth is permitted, the CKS solution 
for straight incident strings evolves such that the bridge 
grows at the constant rate specified by Eqn. ((3j, with the 
bridge orientation 9 and velocity u constant also. Kinks 
travel out along the initial strings, as shown in Fig. [T] 
and the geometry simply scales in size with the time t 
since the collision. 



III. U(l)xU(l) DUAL ABELIAN HIGGS MODEL 

We now explore these dynamics from a field theory 
perspective, including strings which have finite width, 
interact strongly in the region close to the Y-junction and 
may radiate. To do so we employ the dual U(l) model 
of Ref . [30 ] , which has also been used in horizon- volume 
simulations in Ref. [25| • This involves two Abelian Higgs 
models coupled only via the potential term and having 
Lagrangian density: 
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We have followed the notation of Ref. [30] and defined 
the gauge-covariant derivatives as: 
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TABLE I: The energy per unit length of a static string with 
a winding of 2nm in the phase of <j) and of 2irn in the phase 
of ij), with parameters Ai = A2 = 2e 2 = 2g 2 = 2, 77 = v and 
K = O.4VA1A2. Note that an Abelian Higgs string of unit 



winding would yield fi — 2-kt\ 36 



For k — it is well-known that string solutions exist for 
each half of the model [35| (see Refs. 0,0] for reviews). 
These are characterized by the phase of <j> (or tj;) having a 
net winding of 2inn (or 27m) around any closed path that 
encloses the string. For k^0 static and straight (m,n) 
string solutions, involving both halves of the model, were 
found in [301 ]. By way of an example, consider v = r\, 
e = g and Ai = A2 so that /tt(o,i) = Md.ol and also let 
2e 2 = Xi (so that the Bogomol'nyi limit 36] applies in 
both halves of the model). Then if k = OA^/XiX^: 



M(i,i) = °- 840 [M(i,o) +M(o,i)] 



(11) 



and hence parallel /i(o,i) and /i(i,o) strings can reduce 
their energy by combining to give a composite ^(1,1) 
string. Smaller (positive) values of k yield a lower re- 
duction in energy by composite formation, while n > 
0.5 a/A^X^ results in the model being unphysical, because 
then the potential is unbounded from below. Finally, 
negative k yields an increase in energy so composite so- 
lutions are unstable and such values are not of interest 
here. The values of /X( m n ) for the strings involved in our 
simulations are shown in Table [H 
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mretric field strength tensors as: 




i 1 pi/ = (j^A v — o v A^, 


(7) 


Ffxv — d l _ i B v — dyB^. 


(8) 



The only coupling between the two otherwise indepen- 
dent Abelian Higgs models is via the final potential term 
of Eqn. (J4j) . Its form is chosen to ensure that the lo- 
cal U(l) symmetries associated with each Abelian Higgs 
model are preserved: 
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(10) 



since these are directly related to the presence of string 
solutions. 



IV. SIMULATION METHOD 

A. Evolution algorithm 

In order to simulate the collision of two straight strings 
we represent the fields of the model using only their val- 
ues at discrete points in space and time, and then write 
approximations to the second order dynamical equations 
in terms of the fields at these points. There is no unique 
way to proceed, however, for the Abelian Higgs model in 
Minkowksi space-time, it has become popular to use the 
approach of Ref. [33]. In that method, a discrete Hamil- 
tonian is constructed and then the equations of motion 
for the discretized variables are obtained from it. The 
Hamiltonian is chosen to preserve the U(l) gauge symme- 
try of the Abelian Higgs model, at least in a certain dis- 
crete form. As such the simplifying gauge choice Aq = 
may be chosen but the discrete Hamilton's equation for 
Ao, which is the analogue of Gauss' law in the model, is 
preserved precisely by the remaining discrete equations. 



This would be difficult to achieve via the direct replace- 
ment of the derivatives in the dynamical equations with 
finite differences. Note also that this approach may be 
generalized to flat Friedman-Robertson- Walker (FRW) 
cosmologies via the use of a discretized action, as in Ref. 

Since the coupling between the two Abelian Higgs 
models in the present action is via the potential term 
only, and involves no field derivatives, then this does not 
greatly affect the discretization of the system. We there- 
fore applied the above procedure to arrive at our evolu- 
tion algorithm for the U(l)xU(l) model in Minkowksi 
space-time (although our program was derived directly 
from the FRW code of Ref. [38]). Parallel computation 
was made available via the use of the LATfield library 
|39| . with simulations performed across up to 32 proces- 
sors of the UK National Cosmology Supercomputer [40] • 
The (scalar 1 ) fields were represented on the sites of a 
cubic lattice of spacing Aa; and therefore yielded a uni- 
form (as opposed to adaptive) spatial resolution, while a 
constant timestep At was also employed. 



B. Initial conditions 

We desire to start the evolution with two straight, in- 
finite strings moving towards each other, but unfortu- 
nately there are no such analytical solutions known. Us- 
ing the code written for [30]. however, we can rapidly 
obtain numerical solutions for isolated static strings with 
given winding numbers, which we then employ here to 
construct the initial conditions following a procedure sim- 
ilar to that used by Ref. [37| for the Abelian Higgs 
model 2 . 

The isolated string code solves for the radial profile 
of a string using a cylindrical polar coordinate system 
around its centre, which we denote as (r',6',z') and the 
dashes denote that we are in the string rest frame, rather 
than the simulation frame. Via the appropriate choice of 
gauge, the profile can be written as: 
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The gauge field components arc represented half-way along the 
links between the sites (A x on a link parallel to the x-axis, A y 
on a link parallel to the j/-axis, etc.), explicitly transporting the 
phase of the scalar fields across the links. 

See also Ref. [4l| for an alternative approach for the Abelian 
Higgs model, employing instead the Lorentz gauge, and Ref. [42)1 , 
which in fact employs a Lagrangian similar to Eqn. J4}, but in 
a different regime in order to model superconducting strings. 



with the other components of the gauge fields simply 
zero. That is, for a given m and n this code returns the 
four functions /(?"'), a(r'), p(r'), and b(r') for an array 
of discrete r' values. 

In order to obtain the solution for the string in the sim- 
ulation frame, that is a string moving at a speed v in the 
x-direction, we must follow a similar perform a Lorentz 
boost. However, we also wish the string to be rotated 
through an angle a in y-z plane, as shown in Fig. [1] 
We therefore perform a Lorentz boost and simultaneous 
rotation, giving cj> and A^ in the simulation frame as: 



<t> 


= <t>', 
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= <y(A' -vA x ,), 
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= j(A' xl -vA' ), 
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Ay 


= A' y , cos a — A' z , sin a, 


(19) 


A z 


= A y , sin a + A z , cos a. 


(20) 



Throughout this section equations for <fi and A^ will have 
obvious tp and B^ analogues and we will therefore not 
repeatedly note that the second half of the model is to 
be treated in the same manner of the first. From Eqn. 
(fT7)l it can seen that if the above rest frame solution 



with A' n = and r'A' , 



-A'ai sin 9' is inserted, then 



the resulting Aq is non-zero. That is, the solution is not 
of the appropriate form for use in the Aq = evolution 
algorithm. 

Therefore, before applying the boost we perform a 
gauge transform muiA such that A' , becomes zero: 



1 



(21) 



d x i(muJA) = -A' x ,. 
A solution to this equation is given by: 

u, A (x',y') = -- f X A' x ,(X',y',z') dX>, (22) 
m Jo 

and if the static solution above is inserted, this becomes: 



uj A (x',y') = 



o A' 2 + / 



a(yA" 2 + y' 2 ) dX', (23) 



which we evaluate numerically. Note that the integral is 
finite for all x' and y' since a must tend to 1 as r' — ► oo 
in order to give finite energy per unit string length. 

We choose zero as the lower integral limit because then 
loa{x' , y') is an odd function with respect to both x' and 
y' . As a result lo a needs only to be found in one quadrant 
and can then be easily applied to all four: 
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Having determined u>a, then the gauge transform 
yields: 



0' = nf{r')e im( - 6 ' +UA \ 

A' = -(- a (r')cos9' 
j e \r 



(25) 



dy'UA 



with all other components of A' now zero. The y'- 
derivative of loa can be itself written as an integral but 
a simpler and numerically faster procedure is to calcu- 
late 0Ja(x' ,y' ± Sy') for Sy' much less than the y' grid 
spacing and then find the gradient via a (centred) finite 
difference. The fields in the simulation frame are then: 
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where 



x = j(x — vt), 

y' = ycosa + z sin a. 



(26) 
(27) 
(28) 
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However, it is also required to specify the time deriva- 
tives of the fields since the equations of motion are second 
order. The string solution is simply translating at a ve- 
locity v in the x-direction, therefore these derivatives can 
straightforwardly be obtained as: 



d t <f) = —vd x <f> 



(33) 



where d x (j) can be found from 4>(x ± Sx,y,z); and with 
an analogous equation for the gauge field. 

Obviously, the above procedure generates only a sin- 
gle string and so it is required to superpose two such 
solutions, one with positive v and one with negative v, 
but positioned such that they are initially separated by 
a distance much larger than the string width. Since the 
equations of motion are non-linear there is no precise 
means to do this, however for such separations a good 
approximation is simply to sum the single string gauge 
fields A+ and A~ to give the total: 
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Then for the scalar fields: 
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results in a superposition of complex phases and min- 
imal interference of the core of one string due to the 
approximately constant, near vacuum field of the other 
[23, SI Sil • Using these equations, the time derivatives 
must then superpose as: 
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In this way two straight strings are created such that 
they are approaching each other and, given the form of 
the time dependence in Eq. 1311 their centre-lines will 
collide at time t = 0. That is, the initial separation is set 
by the choice of the simulation start time t s tart- 




Cc:(1,1)+(0,-1)-»(1,0) 



Cs:(1,1)+(0,-1)-»(1,2) 







FIG. 3: Schematic illustrations of the three collision types A, 
B, and C considered here, with the difference between the two 
possible connectivities which exist for type C highlighted: Cc 
(cosine) and Cs (sine). 



C. Boundary conditions 

Since the evolution algorithm uses finite differences to 
represent spatial derivatives, then the update of a field at 
a particular site requires the knowledge of neighbouring 
sites. However, at the simulation boundaries the fields 
are not known for all neighbours and a method must be 
chosen to determine a value for them. Here we simply 
employ our initial conditions code (but with t > £ s tart) 
to calculate these unknown field values [4l| . Hence there 
is a halo of sites surrounding the main simulation vol- 
ume whose values we update using the initial conditions 
algorithm after each timestep. 

Note that our boundary conditions will therefore re- 
flect waves travelling along the strings and that the simu- 
lations are reliable only while the boundaries are causally 
unaware of the interaction between the two strings. We 
therefore must chose our lattice size appropriately to 
yield an adequate time to study the collision aftermath 
before artifacts of the boundary conditions reduce the 
reliability of the simulations. 



V. QUALITATIVE RESULTS 

Before revealing our methods for quantitative mea- 
surements of the post-collision environment we must first 
present our basic qualitative results. We limit ourselves 
to model parameters 2 = Ai = A2 = 2e 2 = 2g 2 and 
K = O.4VA1A2 and also choose equal energy scales for 
the two halves of our model 77 = v. As briefly mentioned 
in Sec. |Hj we consider three different sets of initial string 
windings: cases A, B and C, as shown in Fig. [3] We will 
now discuss each one in turn. 






FIG. 5: Energy density slices through a case A simulation, 
showing the initial burst of radiation as the bridge first forms. 
These images are heavily saturated in the bridge core, as is 
required to resolve the small amount of radiation emitted, and 
are taken from an a = 20°, v — 0.2 simulation at t = 21rj — . 
In the uppermost pane, the bridge appears darker than the 
initial strings because it lies in the x — plane, while the 
initial strings are 4.2r; _1 off this plane, and note that the two 
panes share different spatial scales. 



FIG. 4: Isosurfaces of \<j>\ = 0.5r; and \tp\ = 0.5^ from a type 
A collision: (1,0) + (0,1). The upper two frames are from 
t = — 67/ -1 , showing the attraction of the strings towards 
each other and intersection at negative times rather than the 
Nambu-Goto value of t = 0. The lower two frames are from 
the later time of t — +S0rj~ 1 , showing the established bridge 
and Y-junctions. The initial speeds were v = 0.2 and the 
strings made an angle a = 20° to the z-axis. 



A. (1,0) + (0,1) -» (1,1); 



For the collision of a (1,0) string with a (0,1) string 
to plausibly yield a (1,1) string, we have /ii = /12 and 
R = 0.840 and we do indeed find that a composite (1,1) 
region forms for low values of v and a, with a Y-junction 
at either end and having cosine linkage. 

An example with v = 0.2 and a = 20° is shown in Fig. 
2J The basic expectations from the CKS solution are 
indeed seen in this case, albeit with the two strings at- 
tracting each other slightly before the theoretical collision 
time of t = 0, and as can be seen in the upper half of the 
figure, the central region of the upper string is attracted 
down towards the lower one (and vice- versa). The dip 
(or raise) breaks up into waves travelling in opposite di- 



rections 3 , resulting in low amplitude x-displacement pat- 
terns which move along the strings once the Y-junctions 
have formed. There is also a smaller disturbance trav- 
elling along each string with a displacement in the y-z 
plane, created during the parallel re-alignment of the two 
strings at the intersection. In general such waves may be 
more or less pronounced and, for example, as a is de- 
creased then the t < interaction has a greater z-range 
and the strings then undergo many oscillations as the 
strings align parallel to each other. 

At late times, when any such oscillations have subsided 
and the displacement waves are far from the Y-junctions, 
these junctions do not settle down to take on precisely 
the sharp Y-shape as in the Nambu-Goto case. As would 
be expected for interacting strings of finite width, the 
(1,0) and (0,1) strings curve gradually away from the z- 
axis as we move out from the junction, with the radius 
of curvature being a few times the string width. 

The formation of the bridge is accompanied by a burst 
of low intensity radiation, which is shown in Fig. [5] No- 
tice that the energy density in the radiation is tiny com- 



3 In order to understand this, it may be useful to the reader to 
note that x = x+ (s + t) + x (s — t) is the general solution to 
the equation of motion x = <9 2 x/c9s 2 of Nambu-Goto strings and 
that for an initially static string x_|_(s + t) is equal to x_(s — t). 
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FIG. 6: Incident velocities v and angles a which are seen in 
simulations to yield Y-junctions when /xi = ^2 (i.e. 5 = 0) 
and R = 0.840. Here Y denotes the formation of a bridge 
with Y-junctions at either side, points that the strings passed 
through each other and X denotes that the strings became 
locked together to form an X-junction. The CKS prediction 
is that the region beneath the curve could yield Y-junctions. 



pared to that in the string core and that this represents 
the peak of radiation production. The initial burst soon 
ceases and there is then no reliably resolvable radiation 
produced at late times. 

Keeping a fixed and increasing v reveals that, as ex- 
pected from the CKS calculations, eventually the com- 
posite region no longer forms and the strings merely pass 
through each other. The only trace of the interaction is 
then displacement waves, similar to those seen when Y- 
junctions do form. However, the limiting v for composite 
formation is often somewhat below the CKS prediction, 
as shown in Fig. [6] for this symmetric case. Note that we 
need only consider a < 45° due to the symmetry present 
in the initial conditions. This discrepancy between the 
Nambu-Goto and field theory cases is somewhat more ex- 
treme than reported by Ref. [32| for the Abelian Higgs 
model in the type I regime, something that we will return 
to in our conclusions. We will also investigate the transi- 
tion from bridge-forming collisions to non-bridge-forming 
collisions in Sec. IVIII when we discuss quantitative mea- 
surements from the simulations. 

Interestingly at larger a we find that there are still 
two possibilities for the final state of the system. For 
large speeds it remains the case that the strings pass 
through each other, while for lower values the strings 
become locked together and an X-junction forms, as il- 
lustrated in Fig. [7] These are denoted by an X in Fig. [6] 

The Nambu-Goto solution for four strings connected in 
an X-junction is trivial, given the present initial condi- 
tions and the equal tensions /ii — \ii- The X-junction 
itself is simply static by symmetry and the string lo- 




FIG. 7: Isosurfaces of energy density T = 0.5r/ from a colli- 
sion of type A, when a = 40° and v — 0.2, showing the forma- 
tion of an X-junction. Results are shown for time t — 30r; _1 . 



catcd between it and the kinks must then also be sta- 
tionary (since the junction just reflects incident waves). 
However, the Nambu-Goto equations for four connected 
strings cannot yield any constraints because the connec- 
tivity is put in by hand. 

The field simulations show a very similar situation to 
these Nambu-Goto dynamics, albeit for additional mi- 
nor oscillations, as in the Y-j unction cases. There is 
only a small interaction region, but this will have a low- 
ered energy per unit invariant length. The small amount 
of energy liberated by this, and a larger amount of en- 
ergy liberated by the complete retardation of the inci- 
dent strings, must then go into the increased length of 
the string, since there is again little excitation of radia- 
tive modes. However, since the binding is over only a 
short length of string it would be expected that even a 
small perturbation from these very idealized initial con- 
ditions, such as a low amplitude disturbance travelling 
along one of the strings, would easily break up the X- 
junction. The strings would then separate due to their 
tensions and therefore we do not believe that X-junctions 
would be cosmologically important in this model. They 
may, however, be more relevant in non-Abelian models 

m 



B. 



(2,0) + (0,1) -» (2,1): 



A similar situation exists also for the unsymmetric case 
of a (2,0) string colliding with a (0,1) string to yield a 
possible (2,1) composite. Since S — 0.305 it would be 
expected from CKS that, when a bridge forms, it would 
not be static and would not lie parallel to one of the 
coordinate axes but instead traverse in the same direction 
as the heavier initial string and be orientated closer to it. 
This is indeed the result apparent in Fig. [5] although we 
return to this comparison from a quantitative perspective 
in Sec. EEH 

As for case A, there is a burst of radiation as the bridge 
forms, but again the emission appears to be weak and 
limited to the bridge coalescence phase. The distribu- 
tion of Y-junction formation events across the a-v plane, 
shown for this case in Fig. [91 is also similar to case A. No 
Y-junctions are found when CKS solutions are forbidden, 
but the limiting velocity at a given a is again lower than 




FIG. 8: Isosurfaces of energy density T® = 0.5»7 4 from a type 
B collision: (2,0) + (0,l), with a and v as in Fig. H The 
snapshot is for t = 24rj~ 1 . Of the two initial strings, it is the 
(2,0) string which appears thicker and notice that the (2,1) 
bridge is angled toward this string. 
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FIG. 9: Incident velocities v and angles a which are seen in 
simulations to yield Y-junctions when S = 0.305 and R = 
0.840. Here Y denotes the formation of a bridge with Y- 
junctions at either side, points denote that the strings passed 
through each other and X denotes that the strings became 
locked together to form an X-junction. The CKS prediction 
is that the region beneath the curve should yield Y-junctions. 



the Nambu-Goto dynamics would allow. 

Also mirroring case A, large values of a give X-junction 
formation at low v. Even in this unsymmetric case the 
Nambu-Goto solution for a X-junction is trivial and af- 
ter performing a Lorcntz transform to the rest frame of 
the junction, the strings between the junction and the 
kinks would become static. This appears to be be essen- 
tially what is seen in the simulations although with small 
oscillations as before. 



C. (1,1) + (0,-1) -» (1,0) or (1,2)? 

As shown in Fig. [3l the intersection of a (1,1) string 
with a (0,-1) string can yield Y-junctions of two types 
depending upon the connectivity. Related to this is that 
the negative sign in (0,-1) has no absolute meaning with- 



out further qualification, and here we use the minus sign 
to denote that if a = then the two strings carry op- 
posing -0-fluxes. Hence the naive expectation is that, for 
small a, the ip fluxes will annihilate and a (1,0) bridge 
will form, which would yield R — 0.373. We refer to this 
case as Cc, since it has cosine connectivity with respect 
to a. On the other hand, for a close to 90°, the phase 
of ip will wind around the strings in the same direction 
and it might be anticipated that a (1,2) bridge will form, 
giving R = 0.901. We label this case with sine linkage as 
Cs. For this set of initial windings, therefore, the differ- 
ence between the two connectivities has an impact upon 
the value of R and it is not simply cos a — > sin a and 
cos 9 — > sin 9 in Eqn. $5§ for the bridge growth rate. 

Since |5| > R 2 for a (1,0) bridge, then the Nambu- 
Goto prediction is that, as was shown in Fig. [21 bridge 
formation will occur for all a at large initial velocities 
(v > 0.791). On the other hand, this is not true for the 
heavier (1,2) bridge and the prediction is then very much 
like that for cases A and B. There is, however, a small 
region of overlap, where both Cc and Cs allow bridge- 
forming CKS solutions and the Nambu-Goto equations 
make no predictions as to which connectivity will actually 
be given. 

Indeed it is seen in the simulations that cosine bridge 
formation can occur at very high velocities for all values 
of a tested, as shown in Fig. [TUJ An example of Cc 
bridge configuration is then shown in Fig. 1111 Collisions 
yielding sine connectivity appear, as expected, only in 
the bottom right of the a-v plane, while both Cs and Cc 
solutions are seen to populate the region in which CKS 
allows both solution types. There are no bridges seen 
in the region which CKS solutions are not possible, but 
as with cases A and B, the simulations show that the 
Nambu-Goto dynamics do not yield the boundaries be- 
tween allowed and disallowed regions with total accuracy. 

The match between the CKS predictions and the simu- 
lations results is, however, somewhat closer in the present 
case. This is perhaps because, now differing from the 
first two collision types, the incident strings both have a 
finite ^-winding. Conventional cosmic string lore would 
dictate that if two (0,1) strings collided, they would in- 
tercommute and therefore we have good reason to ex- 
pect a significant interaction in this case, rather than 
the two strings passing through each another. Indeed, 
Fig. [T2] shows the centre-lines of the strings detected 
via the winding of the scalar field phases around lattice 
grid squares, and shows that in this type Cs collision, 
the first step toward the formation of a Cs bridge is an 
intercommutation of ^-strings. The strings in the inter- 
action region are then two (0,1) strings with a (1,0) string 
between them and the system has two paths via which 
to form Y-junctions: either (i) the i/'-strings upzip from 
the binding </>-string and leave a type Cc bridge, or (ii) 
the tp strings zip-up further along the <p string producing 
a Cs bridge. The illustrated case is one in which both 
are possible under the Nambu-Goto approximation, but 
the field theory chooses the latter. Note that there is a 
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FIG. 10: Incident velocities v and angles a which are seen in 
simulations to yield Y-junctions for case C: (l,l) + (0,-l), and 
compared to the CKS predictions from Fig. [2] The final state 
of the system is denoted by a c when a (1,0) cosine bridge and 
Y-junction pair formed, an s when a (1,2) sine bridge was seen, 
an X when a X-junction formed and a dot when the strings 
passed through each other. Additionally (c) denotes a case 
when a cosine bridge formed and grew, but that a displace- 
ment wave caused a late-time intercommutation between the 
(1,1) and (0,-1) strings, after which the bridge collapsed. 



FIG. 12: The reconstructed string centre-lines for a type Cs 
collision: (l,l) + (0,-l)— >(1,0), occurring in the region of the a- 
v plane where the Nambu-Goto dynamics suggests that both 
Cc and Cs are possible: a = 70° and v = 0.2. Results are 
shown for t = — 2477~ 1 (left), t — (centre) and t = 24r? _1 
(right), showing the <f> string centre-line in blue and the %j) 
centre-line in green. Centre-lines are detected using the wind- 
ing of the scalar field phases around lattice plackettes (see Sec. 
IVII and Appendix |A1 . Lengths are shown in lattice units: 
Ax = 0.5r? _1 . 




FIG. 11: Isosurfaces of energy density Tq = 0.5t? 4 from the 
collision of (1,1) string with a (0,-1) string with initial speed 
v — 0.6 and a = 20°, showing a (1,0) bridge has formed. The 
snapshot is shown for t = 4477 -1 . 



further -0-string intercommutation event before the end 
of the simulation but the (1,2) bridge is stable and such 
intercommutation events are inconsequential. 

Even when neither Cc nor Cs bridges form, there is still 
the intercommutation between the '0-strings. However, 
again a second intercommutation event occurs and, now 
in contrast with the illustrated case, the binding energy 
of the strings is not sufficient to hold them together so 



the (1,1) and (0,-1) strings simply separate. 

As in collisions of type A and B, there is a burst of 
low intensity radiation as the bridge forms, but with the 
emission seeming to decay away at late times. 



VI. MEASUREMENT OF THE BRIDGE 
LENGTH, VELOCITY AND ORIENTATION 



A. Symmetric case: hi 



fJ2 



Having seen that Y-junctions form in the present 
model in a very similar manner to the CKS solutions, 
we desire to test the CKS predictions quantitatively by 
measuring the bridge length as a function of time, as 
well as its velocity and orientation. We start by detailing 
our method for case A since then symmetry dictates that 
the only possible (mean) bridge velocity is zero and the 
only possible (mean) orientation is parallel to either the 
z-axis or the y-axis. We restrict ourselves to a < 7r/4 as 
before, thanks to the symmetry present in the problem, 
and from our previous observations need only consider 
bridges forming along the z-axis. 

In order to measure the bridge length in this symmetric 
case we may merely count the number of sites along the 
line x = y = which have \4>\/ri < e or {ipl/v < e, with 
e some chosen threshold. Multiplication by Ax and divi- 
sion by 2 then yields a measure of the physical half -bridge 
length, which for the static bridge is equal to the invari- 
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ant half-length S3 . Note that since no other string lies on 
this line and any initial oscillations along the bridge are 
only transiently important, then this is a reliable method 
which gives easily interpreted results. 



B. General case: ^ 7^ \xi 

Unfortunately the above procedure is not readily appli- 
cable to the general case and we must therefore employ a 
more complex method to determine the bridge length for 
cases B and C, as well as additionally measuring the non- 
trivial bridge velocity and orientation which they yield. 

In past work [15|, [2l|, [3a, H] , Abelian Higgs strings have 
been detected using a net winding of the scalar phase 
around the smallest closed loops resolved in the simula- 
tion: the lattice plackettes. However extending this to 
detect, for example, (1,1) strings is non-trivial since it 
is not guaranteed that the phase of <j> has a net winding 
around precisely the same plackette as the phase of tp. 
Furthermore, for (2,1) strings it is not even assured that 
<f> will exhibit a 47r winding around a plackette rather 
than there being merely two close-by 2tt windings. 

Our approach is therefore to firstly trace the path of 
every (1,0) and (0,1) elemental string through the simu- 
lation volume, such that a 4-7T winding in cj> denotes that 
two (1,0) strings thread the plackette 4 . We then detect, 
say a (2,1) bridge, by searching through these paths for 
regions in which precisely two 4> string paths and one tp 
path string approach within a certain distance of each 
other. 

The easiest method of reconstructing the string paths 
is to suppose that a straight segment of string of length 
Ax passes through each plackette which has a net wind- 
ing. Unfortunately constructing the paths out of an ar- 
ray of perpendicular segments can yield to a significant 
over-estimate in the total length. For the cosmological 
simulations of Refs. [15| and [2l|, this overestimation 
was countered by smoothing the paths on the scale of 
the string width, although here we may employ a quite 
different approach. 

A basic observation from the simulations is that any 
bridge that forms is straight and, even if oscillations are 
present initially, these quickly decay. We therefore take 
the change in position vector Ar along each segment and 
perform the vectorial sum over all detected bridge seg- 
ments. We then measure the physical bridge length 2Z3 
as the modulus of the resulting vector, divided by the 
total winding (|ro| + \n\) of the bridge: 



This effectively fits a straight line through the \m\ + \n\ 
element string paths in the selected bridge region and 
then determines its length 5 (although note that it is in 
fact only a function of where the elemental paths cross 
the bridge detection thresholds). 

However, we additionally perform an interpolation in 
order to more accurately locate the intersection of the 
string centre-lines within the plackettes, for which we em- 
ploy the method described in Appendix [A"l The vectors 
Ar are then taken to link two such intersections. Actu- 
ally for our vector-based ^3 estimator this is important 
only for the ends of the vector and has a small effect. 
However, this interpolation is more important for the 
measurement of bridge velocity u and the angle 9, the 
methods for which are described momentarily. Further- 
more, the use of our interpolation scheme and a direct 
J2 |Ar| sum of lengths yields a second I3 estimator that 
we have employed as a check, and which we find performs 
surprisingly well. 

In principle, the angle 9 could be taken as that between 
the z-axis and the summed bridge vector J2 Ar, however 
this would potentially yield a small systematic error. For 
example, where a (2,1) bridge splits into (2,0) and (0,1) 
strings, the vector will be biased towards the (2,0) string 
and, as a result, 9 will be over-estimated. Although the 
effect will be small, the CKS predictions for 9 show a 
very mild dependence upon v for fixed a, R and 5, and 
we wish to accurately explore this in the field theory case. 
We therefore perform a second vectorial sum using only 
the central 80% of the bridge and then measure 9 using: 



tan 9 



Sgo% Ar v 

1^80% ^ r z 



(39) 



that is, the angle between this new resultant vector and 
the z-axis. 

Finally, in order to determine the bridge x-velocity u, 
we use the weighted mean of the segment cc-coordinates 
over the central 80% of the bridge: 



T,80%( X ++ X -)\ Ar \ 

2E|Ar| 



(40) 



2/3 



1 



l£ Ar l 



(38) 



where x± are the endpoints of a segment. Note that the 
sub-plackette interpolation scheme is more important in 
this case since it derives from the precise values of x± 
even in the central regions of the bridge, rather than just 
at the extremes. 

Given these measurements of ^3 and x, we determine 
dls/dt and u by performing linear fits for late times, while 
for 9 we take the mean during the corresponding period. 
Note that we do not combine our measurements of u and 
dl^/dt to yield S3, but rather compare our measurements 



4 Note that we use the gauge-invariant winding measure of Ref. 
[44J since in the discrete case the winding in the scalar phase can 
be removed by a finite gauge-transformation. 



Both rj> and i/i-strings are taken to traverse the bridge in the 
same direction and therefore a (1,1) string, for example, would 
not yield erroneous cancellations in the vectorial sum. 
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of physical length to the CKS predictions in this quantity 
using Eqns. flU and ©. 



VII. NUMERICAL RESULTS 

We now present the results of applying the above al- 
gorithms to the three cases for which we presented qual- 
itative results in Sec. [Vj As before, our simulations all 
have 2 = Ai = A 2 = 2e 2 = 2g 2 , rj = v and k = 0.4-\ZAiA 2 . 



A. (1,0) + (0,1) -» (1,1) 

In the symmetric case of a (1,0) string colliding with 
a (0,1) string, we find that the bridge half-length varies 
as shown in Fig. [13] for a = 20° and v = 0.2. The 
CKS prediction is also indicated on the plot, being the 
straight line for which I3 = when t = 0, and the two sets 
of measurements from the simulation do approximately 
track this. The difference between the latter two is simply 
due to the choice of a different threshold e below which 
both |<^>|/77 and \i>\jv must be in order for a site to be 
considered as part of the bridge. Of course, it would be 
expected that the case of e = 0.75 would show larger 
I3 values than e = 0.5 since a larger threshold will be 
crossed further from the string centre-lines. 

If the CKS solution was precisely followed by these 
centre-lines, it would be expected that these two mea- 
sures each would show I3 = ts^ KS + c, where c is con- 
stant and equal to the distance between the point on the 
z-axis at which the threshold is crossed and the point 
where the three string centre-lines meet. This is approxi- 
mately what is seen, although it should be noted that the 
collision of the centre-lines occurs slightly before t = 
due to the attraction between the strings. Of course, for 
t >> r/ _1 the offset c will become negligible and hence 
from a cosmological perspective we are really only in- 
terested in whether the late-time gradient is accurately 
predicted by the CKS solution. 

The measured gradients for various a and v values 
are shown in Fig. Q3] and compared to the CKS pre- 
dictions. An approximate uncertainty estimation is per- 
formed such that dls/dt is taken from a linear fit to the 
final third of the apparent linear region, with the error 
bar shown being the standard deviation across the three 
thirds. This method is sensitive to both systematic dif- 
ferences between the early- and late-time dynamics and 
to the measurement uncertainties. The plot shows that 
when Y-junctions form there is excellent agreement be- 
tween the simulations and the Nambu-Goto predictions. 
However if v is increased at fixed a, then there is a cer- 
tain critical value v c where dl^/dt drops suddenly to zero 
and away from the CKS solution. There is not, as one 
might have expected, a gradual divergence from the CKS 
predictions, and the results seen in Fig. [5] are not due to 
dls/dt in our model slowly falling away from the Nambu- 
Goto value and reaching dl^/dt = at a lower value of 




FIG. 13: The measured Z3 as a function of t measured from 
a Aa; = 0.57? -1 simulation for a case A collision: (1,0) + (0,1), 
with a — 20° and v = 0.2. Results shown are derived from 
the count of sites on x — y = with \cj>\ < erj and \ip\ < eu and 
using e = 0.5 (blue crosses) and e = 0.75 (green triangles), 
while the CKS prediction is shown by a dashed grey line. Blue 
cicles indicate the results for e = 0.5 from a shorter simula- 
tion but with Ax halved from 0.5?? -1 to 0.25»7 , highlighting 
that the simulated dynamics are not precisely those of the 
continuum, while the measurements themselves are accurate 
to within Ai/2 and the corresponding uncertainties are too 
small to be shown on this plot. 




FIG. 14: The measured dl^/dt values from a class A collision, 
shown as a function of v for fixed a, compared to the CKS 
predictions. Results are shown for a — 10° (uppermost), 20° 
(middle), and 25° (lower). The simulations had Aa; = 0.5 and 
were of sufficient size that signals emitted from the box centre 
at t — would reach the corners of the yz-plane by t — 60r; _1 , 
except for (a = 20°, v = 0.2), (a = 20°, v = 0.3405) and 
(a = 25°, v = 0.31) which ran till t = 120, 90 and I2O77- 1 . 
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v. Instead it is the case that: either Y-junctions do not 
form at all, or they form and an dl^/dt value equal to (or 
at least extremely close to) the CKS value is observed. 



B. (2,0) + (0,l) -» (2,1) 

The measurements of dl^/dt, u and 9 from a collision 
of type B are shown as a function of v for fixed a = 20° 
in Fig. 1151 with uncertainties estimated via the same ap- 
proximate method as for case A. Results indicate that, 
as in the type A collisions, if Y-junctions form then these 
three measures are well-predicted by the Nambu-Goto 
physics. Even the slight departure from u ex v, a case 
which would appear on our u/v plot as a horizontal line, 
is quite obvious and in concordance. However, resolving 
the change in 9 with v, at fixed a, is more challeng- 
ing, being just 0.6° across the entire range over which Y- 
junctions are seen to form. Hence our simulations have 
struggled to concretely resolve the trend, although they 
do give an indication of the slight increase at larger v. It 
did not seem worthwhile to utilize very large simulations 
in order to reduce the measurement uncertainties since it 
is surely of little consequence whether the CKS predic- 
tions get the bridge orientation wrong by small fractions 
of a degree, and in any case, we explore this variation 
more completely for case C. 

As in case A, the plot shows that when v is increased to 
a certain critical value, the CKS predictions suddenly fail 
and Y-junctions do not form for greater speeds. Again 
there is no gradual reduction in the late-time dl^/dt 
value. However in this case, with symmetry not fixing u 
and 9, we see additionally that these two appear to match 
the CKS predictions right up to this critical speed. 



C. (l,l) + (0,-l) -» (1,0) 

Case Cc provides a quite different collision, in that 
1 5 1 > R 2 , and also the initial strings both have a net 
winding in ip. As shown in Fig. 1161 the dl^/dt results 
are for type Cc collisions at a = 60° are in complete 
concordance with the CKS predictions. This agreement 
appears to be maintained also for both u and 9, with the 
more sensitive dependence of the latter upon v, relative 
to case B, clearly resolved by these data. 



VIII. DISCUSSION OF THE RESULTS 




FIG. 15: The measured physical bridge half-length growth 
rate dls/dt, bridge speed u and bridge orientation 8 from type 
B collisions with a = 20°, compared to CKS predictions. 
The simulation had Aa; = 0.5 and were of sufficient size that 
signals emitted from the box centre at t — would reach the 
corners of the j/z-plane by t = 60?7 _1 (except for v — 0.39 
when this value was 120?? -1 ). 




FIG. 16: The measured physical bridge half-length growth 
rate dl^/dt, bridge speed u (blue) and bridge orientation 6 
(red) values from type Cc collisions with a = 60°, com- 
pared to CKS predictions. The simulations at low speed had 
Aa; = 0.5, which was then decreased to accommodate Lorentz 
contraction, with the simulation size such that signals emit- 
ted from the box centre at t = would reach the corners of 
the yz-plane by t = 60?? _1 . 



As noted earlier, the Nambu-Goto action is not valid 
in the vicinity of a Y-junction and ignores, for exam- 
ple, the attraction between the strings in that region. In 
this section, however, we present a discussion of why we 
would in fact expect the Nambu-Goto dynamics to be a 
good description of the late-time behaviour seen in our 
simulations. To do so we first re-derive the CKS solution 



using energy and momentum conservation for Nambu- 
Goto strings, rather than the equations of motion, and 
then we discuss the changes that the field theory case 
would yield. A corollary of this will be a direct physi- 
cal interpretation of the CKS constraints on Y-junction 
formation from the intersection of strings. 
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FIG. 17: The geometry involved in the tension-based inter- 
pretation ol Eqn. (J5J. The Y-junction is at position Y and 
the two kinks are at Kl and K2, while the region upon which 
the tensions in question act is highlighted in gray. The dashed 
blue lines indicate the consumption of the string 1 due to the 
motion of Y. (Note that both the unit vector d along the 
line YK1 and the orthogonal velocity w, have in general fi- 
nite components in the x-direction, therefore the projections 
of them onto the y-z plane would not normally be perpendic- 
ular, although we shown them as orthogonal for illustrative 
purposes) 



Energy-momentum derivation of the CKS 
solution 



In order to derive the CKS solution from energy 
and momentum conservation, we must first consider the 
energy-momentum tensor T£ of a Nambu-Goto string 
with energy per invariant unit length fi. Labelling space- 
time coordinates as (t,x,y,z) and considering the string 
to lie along the z-axis while travelling at speed v in the 
x-direction, then: 



T£ = M 



7 


—vj 





J 7 


w 2 7 




















l/ 7 



5(x)6(y). (41) 



That is, the x-momentum per unit physical length is (ivy 
and so per unit invariant it is just (iv. The term T* is the 
tension of the string, and it is important to note that this 
is reduced by a factor 7 for a moving string and is not in 
general equal to the mass per unit invariant length. 

We will allow some guiding input from the CKS solu- 
tion and assume that the system has the geometry shown 
in Fig. Q~71 with straight strings except for kinks at Kl 
and K2. We therefore write the equation for the line YK1 
as: 



w£ 



(42) 



where w is the velocity of the string, d is a unit vector 
pointing from Y to Kl and a is the invariant length mea- 
sured from the point wt. Note that we take w ■ d = 0, 
which is effectively a choice of gauge. 
The Y-junction lies at position: 



ry = wt • 



CTy 



U 

s 3 sin 0/tu 
s 3 cos 0/7„ 



(43) 



If Y moves such that cry becomes positive, then invariant 
length is removed from string 1. Following the notation 
of CKS we label the production rate of string 1 at the Y- 
j unction as si and hence cry = —Sit. We also have that 
Kl travels at the speed of light and therefore its position 



i"Ki =wt-\ d = I sin a/7 I t. 

> cos a/7 



(44) 



The value of ctki is given by the invariant string length 
traversed by the kink per unit time and is simply ctki = t. 
Therefore the invariant length of the line YK1 is just 
(<tki - cry) = (1 + h)t- 

In principle, from these line end-points, we can now 
find expressions for s\, w and d in terms of S3, u, and 
9; while applying a corresponding procedure for the line 
YK2. However, for the present discussion we only need 
to determine the total z-momentum from these two lines. 
For the first of them this is revealed by the addition of 
the last two equations as /hi (|43"f +^iisi (|44 |) . resulting in: 

YKl / • s • cos# . cos a , . 

P z = (ii(l + si)tw z = (iis 3 t + (j,isi 1 (45) 

7u 7 

Hence the total z-momentum from both YKl and YK2 

is: 

S3 cos 9 . . . .cos a . . 

p 2 = (Mi+/i 2 ) 1 + (fiisi + (i 2 S2) 1. (46) 

7u 7 

This can be simplified by noting that conservation of en- 
ergy at Y implies that: 



M1S1 +M2S2 +M3S3 = 



(47) 



in the Nambu-Goto case, since there is no radiative emis- 
sion. Therefore we can write the combination of si and s 2 
in terms of S3, yielding the rate of change of z-momentum 
as: 



cos 9 cos a 

Pz = s 3 (Mi + M2)— M3 



7« 



7 



(48) 



That is, the growth of the bridge necessarily leads to an 
accumulation of z-momentum along these two lines. 

Since the incident strings in the regions beyond the 
kinks simply continue with velocities ±w in the in- 
direction, and symmetry requires the bridge to travel 
in the x-direction also, then the only z-momentum in 
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the highlighted region of Fig. [17] is that due to the 
strings along YK1 and YK2. This is provided by the 
^-components of the tensions that act externally on the 
highlighted region: 



cos a cos 9 

F z = [Hi + H2)—, M3" 



7 



7« 



(49) 



and therefore we immediately have a physical condition 
for bridge growth, simply: 



F z >0. 



(50) 



Note that this is not the sum of tensions at the Y-junction 
itself, which might naively be expected to be the relevant 
tensions for the growth condition. It is also not necessar- 
ily true for the tension components along the direction 
in which bridge growth actually occurs. 

Of course, if we now set F z equal to the rate of change 
of momentum, then we trivially obtain Eqn. ([3]). That 
is, we may derive the expression for S3 based upon the as- 
sumptions that the string energy and momentum are con- 
served, straight strings join the Y-junction to the regions 
causally disconnected from it, the bridge is straight 6 , and 
that the energy-momentum tensor is that of Eqn. (|4i"j) . 



B. Explanation of the accuracy of the Nambu-Goto 
approximation 

The above assumptions will not be precisely met in 
the field theory case, but now suppose that, even then, 
the geometry of the system at a certain time t is approx- 
imately as shown in Fig. [TTJ However, allow for dif- 
ferences relative to the Nambu-Goto case including that 
the kinks at Kl and K2 are smoothed on a scale close 
to the string width and that also there are displacement 
waves in their vicinity, which were left over from the ini- 
tial bridge formation. Suppose additionally that there 
will be a significant attraction between the strings in the 
region close to the Y-junction and that their paths will 
gradually curve towards each other rather than there be- 
ing a sharp Y-shape. 

So long as t >> ?/ _1 , these changes will have little ef- 
fect on the large-scale geometry. Therefore we may start 
to follow the above discussion. We must, however, now 
ask ourselves whether energy-momentum conservation at 
Y would take the same form as Eqn. ([47]) . Even with the 
interaction region close to Y, if the shape of the strings re- 
mains constant in time and simply translates with speed 
dls/dt, then the energy associated with the interaction 
would not change. Additionally, if there is merely trans- 
lation of this region, there is no excitation of radiation. 



Note that if the bridge is straight for all time then it must have 
a velocity that is uniform along its length and the velocity must 
then be in the x-direction by symmetry, hence this is not an 
independent assumption. 



Hence it appears that the energy conservation equation 
could be unchanged. 

Therefore, while the z-momentum in the shaded region 
of Fig. [17] would not follow precisely Eqn. (|48ll , for late 
times this difference might be negligible. Further, the 
tensions external to the shaded region would be simply 
the same as in the Nambu-Goto case (for a particular 
and u) and hence equating their z-componcnts with the 
rate of change of z-momentum, we then would obtain S3 
as being very close to Eqn. ([3]). Therefore we expect 
that a plausible late-time solution for the field theory 
case is one which tends towards a CKS solution, despite 
the interaction near Y. That is, of course, assuming that 
a physical CKS solution with S3 > exists for the given 
a and v. 

However, these arguments do not guarantee that the 
Nambu-Goto dynamics will be a good description of the 
system and the initial distortion of the strings due to 
their attraction and the associated oscillations mean that 
the strings do not initially follow the above geometry all 
that closely. There is also an initial burst of radiation, 
as well as the finite width of the strings, both of which 
complicate the field theoretic case. As we have seen, the 
CKS solution appears to be only a reliable indicator of 
when Y-junction formation is not possible, and of the 
late-time dynamics when it is. 



IX. CONCLUSIONS 

Our results indicate that for the present field theory, 
involving two coupled Abelian Higgs models, the CKS 
solutions for Nambu-Goto strings give largely accurate 
predictions as to whether or not Y-junctions will form 
when two straight strings collide. Like the results found 
in Ref. [32| for the Abelian Higgs model in the type I 
regime, we find that they are not entirely accurate, how- 
ever, for two of the three initial string pairings studied 
here we see a noticeably poorer match than in that ref- 
erence. We believe that this difference is due to the fact 
that for our third collision type, both incident strings had 
finite ^-winding and therefore an initial intercommuta- 
tion interaction was highly likely, even with the presence 
of the 0-half of the model. This initial interaction ap- 
pears to aid the alignment of the strings and the initial 
formation of the Y-junction. This observation is further 
compounded by noting that in the Abelian Higgs model 
studied by Ref. [32|], the bridge formation is also pre- 
ceded by an intercommutation event. Hence, we note 
that the applicability of the CKS bridge formation pre- 
dictions to field theory models is likely to be sensitive to 
the model employed, and as we have seen, to the exact 
nature of the strings involved. 

However, in all three of our collision types we find that 
when Y-junctions did form, the late-time dynamics of 
the system was very accurately described by the corre- 
sponding CKS solution. Further, given the discussion 
in the previous section, we believe that the model de- 
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pendence of the initial bridge formation is likely to be 
less significant for the late-time dynamics. Therefore, 
despite the break-down of the Nambu-Goto approxima- 
tion near the Y-junction, the CKS approach appears to 
be a very powerful method for studying Y-junctions in 
local string models. On the other hand, global strings, 
with significant long range interactions are likely to be 
poorly represented by the Nambu-Goto dynamics. 

We find that, in our model, bridge formation is ac- 
companied by a short period of weak radiative emission, 
but find little evidence of such emission at late times. 
Hence this suggests that radiative decay due to bridge 
creation is unlikely to be an effective means of dissipat- 
ing the energy in a cosmological string network. Our 
results therefore offer little comfort to authors who fear 
that networks of string capable of bridge formation might 
"freeze out" and, in contradiction with observation, grow 
to dominate the universe. While our radiative emission 
conclusions are only relevant for the model studied here, 
our arguments from the previous section also suggest that 
local strings in general will mirror the results here. 

A Nambu-Goto simulation using the CKS method 
might well assist the simulation of bridge-forming mod- 
els over horizon scales, since in principle it would avail 
a greater dynamic range than is possible in field theory 
simulations. This is more relevant for these superstring- 
inspired models than for traditional U(l) cosmic strings, 
since the computational outlay for the additional fields 
present is greater, but also the complex network dynam- 
ics require study over a longer time period. However 
with any such simulations, great care must be taken to 
note the differences in the results between Nambu-Goto 
simulations for traditional strings and the U(l) field the- 
ory counterparts, which automatically include a greater 
depth of physics [15|, [20, [2l[ . These differences are, how- 
ever, largely at small scales and do not preclude the 
usefulness of Nambu-Goto simulations with Y-junctions. 
The long term aim must be, of course, to link the string 
network properties to observations and to assess the dif- 
ference between the signatures of traditional U(l) strings 
and cosmic superstrings. 
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APPENDIX A: SIMPLE SUB-PLACKETTE 
INTERPOLATION METHOD 

Given a net winding of the phase of <j> around a lattice 
plackette, there is no unique means of using the mag- 
nitudes and phases of the corresponding scalar field at 
the four plackette corners in order to determine to loca- 
tion of </> = on the plackette surface. In principle one 
may use the known string profile to perform a best-fit de- 
termination of the position and orientation of the string 
centre-line as it intersects the plackette, but for present 
purposes we do not require such a complex procedure. 
Here we use a simple and computationally rapid method 
that, importantly, is guaranteed to yield = coordi- 
nates that lie within the plackette grid square. That is, 
we estimate the coordinates as: 



E l= i x *l 



-l/m 



Ei=l \<Pi 



-1/t. 



(Al) 



where x» are the positions of the four plackette corners, 
(pi the field values at them, and 2-Km is the net winding 
around the plackette. The index —l/m is included in 
this expression because, close to its axis, an ideal static 
(to, n) string has: 



|&|ocr" 



(A2) 



where r is the radial coordinate. This is exactly as in the 
Abelian Higgs model, even though there is the additional 
coupling in the present case. The x expression above 
therefore involves an approximate r _1 weighting of the 
plackette corners so that x will be drawn towards to those 
which appear to lie closer to the string centre-line. 

A clear downside to this method is that, for a plackette 
with sides in the x and y directions, the value of x varies 
with the true y-coordinate of the = point, however, 
it is essentially a zero-cost means of improving the string 
paths through the simulation volume. Unlike a smooth- 
ing operation, this process encapsulates greater informa- 
tion from the fields in the simulation and, for the example 
of a straight string of constant ^-coordinate, improves our 
knowledge of that cc-coordinate, which smoothing cannot 
achieve. It is certainly sufficient for its application in the 
present article. 

A by-eye assessment of its usefulness is afforded by Fig. 
[T2l in which a series of 90° (or 45°) steps is not seen and 
instead the centre-line paths are relatively smooth. It is 
of course the case, however, that some obtuse kinks are 
visible. 
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